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The hypermetric cone is denned as the cone of semimetrics 
satisfying the hypermetric inequalities. Every Delaunay poly- 
tope corresponds to a ray of this polyhedral cone. The De- 
launay polytopes, which correspond to extreme rays are called 
extreme. 

We use this polyhedral cone and the closest vector problem 
to present a new technique that allow to find, from a given 
extreme Delaunay polytope, some new ones. 

Then, we show some examples of applications of this tech- 
nique in low-dimensions. 

A distance vector (dij)o<i<j<n € R N with N = is called an (n + 1)- 

hypermetric if it satisfies the following hypermetric inequalities: 

n 

H{b)d = hbjdij < for any b = (6i) <i<n G Z n+1 with h = 1 . 

0<i<j<n i=0 

The set of distance vectors, satisfying all hypermetric inequalities, is called 
the hypermetric cone and denoted by HYP n+ \. 

In fact, HYP n+ \ is a polyhedral cone (see (7) p. 199). Lovasz (see (7) 

p. 201-205) gave another proof of it and bound max \h\ <n!2"( 2 n n ) for any 
vector b = (6j)o<«<n-i defining a facet of HYP n . 

There is a many-to-many correspondence between Delaunay polytopes with 
p vertices and the elements of the hypermetric cone HYP p (see (7)). So, the 
study of Delaunay polytopes is equivalent, in a way, to the study of hyperme- 
trics. The Delaunay polytopes, whose corresponding face is an extreme ray, 
are called extreme. 

Two extreme ray of a polyhedral cone are called adjacent if they generate 
a two-dimensional face of this cone. Therefore, it induces a natural adjacency 
relation between extreme Delaunay polytopes. Our main purpose here is to 
explain how, given an extreme Delaunay polytope, one can find the adjacent 
extreme Delaunay polytopes. 

This technique allow us to find some new extreme Delaunay polytopes and 
is illustrated in some low-dimensional cases. 

The following closest vector problem is a classic problem of discrete geometry 
and is heavily used in our computations: 
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Closest Vector Problem: 

(i) Given a lattice L C R n , a vector x £ R n and r > 0, test if there exist 
v £ L such that \\v — x\\ < r. 

(ii) Given a lattice L C R n , a vector x € R n and r > 0, find all the vectors 
v G L such that \\v — x\\ = r. 

This classic problem is NP-hard (see (8)). There are several software pack- 
ages doing this computation (see (11) and (14)). 

1. The hypermetric cone 

For more details on the material of this Section see Chapters 13-16 of (7). 
For the use of afhne basis and hypermetrics for computing combinatorial types 
of Delaunay polytopes, see (9). 

Denote by S(c, r) the sphere of center c and radius r. For every family 
A = {^o, • • • , v m } of vertices of a Delaunay polytope P circumscribed by the 
sphere S(c,r) one can define a distance function d& by (d A )ij = H^i — v j\\ 2 - 
The function d& turns out to be a hypermetric by the following formula (see 
(1) and (7) p. 195) : 

m 

H{b)d A = btbjidA)^ = 2(r 2 - \\Y,b t v t - c\\ 2 ) < . 

0<i,j<m i=0 

On the other hand, Assouad has shown in (1) that every d € HYP n+ \ can be 
expressed as g?a with A being a family of vertices of a Delaunay polytope P 
of dimension less or equal to n. 

Definition 1. Let P be an n-dimensional Delaunay polytope with vertex-set 
V. 

(i) A family vq, ■ ■ ■ , v n of vertices of P is called an affine basis if for all 
v € V there exists an unique vector b = (6j)o<i< n G Z n+l , such that 

n n 

biVi = v and b% = 1 . 

i=0 i=0 

(ii) The Delaunay polytope P is called basic if it has at least one affine 
basis. The vertices of an affine basis are called basic vertices. 

All Delaunay polytopes of dimension less, than 6, are basic (see (9)). The 
adjacency method, described later, applies only to basic Delaunay polytopes. 
There is no proof that all Delaunay polytopes are basic and the only indication 
for this conjecture is that all known Delaunay polytopes are basic (see (7)). 

The metric cone MET n+ \ is defined as the set of functions d = (dij)o<i<j<n 
satisfying all triangular inequalities dij < dik + djk', such inequalities corre- 
spond to the hypermetric inequality with vector b, such that 6j = bj = 1, 
bk = — 1 and bi = 0, otherwise. Therefore HYP n+ \ C MET n+ ±. 
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Given d £ HYP n+ \ we define 

n 

Ann(d) = {b G : ^ h = 1 and = 0}. 

i=0 

Theorem 1. Let P be an n-dimensional extreme Delaunay polytope and B = 
{vq, . . . ,v n } an affine basis of P. Then the following statements hold: 

(i) The mapping b — > J2?=o^i v i establishes a bisection between Ann(d) and 
the vertex-set of P. 

(ii) It is possible to compute (in finite time) the list of hypermetric inequal- 
ities H(bi)d < with 1 < i < m incident to d&. 

(in) P has at least ("-+ 2 K n+1 ) _ \ vertices. If it has exactly ("+ 2 K n + 1 ) _ \ 
vertices, then all bi correspond to facets of HYP n+ i . If P has more vertices, 
then it is possible to compute in finite time which of the bi correspond to facets 
ofHYP n+1 . 

Proof. Let S(c, r) be the empty sphere around P and let us write the affine 
basis as t>o, . . . , v n . The equality H(b)d& = corresponds to 

H(b)d B = 2(r 2 -\\^ =0 b i v i -c\\ 2 ) = 

= 2(r 2 - \\v - c + ELi bi(vi - v )\\ 2 ) = . 

This establishes a bijection between solutions of the equation H(b)d& = and 
vertices of the Delaunay polytope P. 

Moreover, the solution-set is exactly a closest vector problem; therefore, it 
is solvable in finite time. Note that all = (0, . . . , 1, . . . , 0) with < i < n 
belong to Ann{d&); for those vector H{ei) = 0, therefore we do not list them 
in the list of vector (6i)i<i< m . 

The Delaunay polytope is extreme; therefore, the rank of the system of 
equations H(bi)d = is equal to "^" 2 +1 ^ — 1. This implies m > "^" 2 +1 ^ — 1. 
Since the vertex-set contains also the affine basis, it has at least n + 1 + m = 

(n+l)(n+2) _ ^ ver ^j ceg _ 

Any extreme ray of HYP n+ i is incident to at least n ^ n ^ 1 > _ i facets. There- 
fore, if P has ( n+1 K"+ 2 ) _ i vertices, then all the H{bi) are facets of HYP n+ i. 

Let us assume that a facet H(bi) is redundant. Then, H(bi) can be expressed 
as a sum with positive coefficients of other facets of HYP n+ \: 

H{bi) = ^ j a j H(b' j ) with aj > 

The relation H(bi)d-B = implies that for all j E J, we have H(b'^)d^ = 0. 
Therefore, b'j belongs to the list of bi. So, the problem of redundancy is reduced 
to the one for the list of facets H(bi) with 1 < i < m, which is solvable by 
linear programming (see polyhedral FAQ 2 ). 

□ 

* http:/ /www. ifor. math. ethz.ch/~f Likud a /polyfaq/ polyfaq.html 
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Above theorem is very useful for finding hypermetric facets. The Delaunay 
polytopes, whose number of vertices is exactly ( ra+1 K"+ 2 ) _ \ are especially 
useful for the adjacency method, by the simplicity of the computation of their 
adjacencies. The Schlafli polytope and the infinite serie of extreme Delaunay 
polytopes in (12) satisfy this lower bound. 

Theorem 2. HYP$ has at least 298592 facets and at least 86 orbits of facets. 

Proof. We know that all facets H(b) of HYP n yield facets H(b') of HYP n+1 , 
with V obtained by adding a zero to b. Since HYPj has 14 orbits of facets 
(see (10)), this yields 14 orbits. 

Moreover, all hypermetric facets of CUT n are facets of HYP n . In (4), a 
conjectural list of orbits of facets of CUT$ is computed. This yields 16 orbits 
of facets of HYP$. 

We know two extreme Delaunay polytopes in dimension 7: the Gosset poly- 
tope and a polytope found by Erdahl and Rybnikov (see (12)). By computing 
all affine basis of them and using above theorem we obtain some more facets. 
It is known (see (7) p. 229) that the switchting of a facet of HYP n is again 
a facet of HYP n . Combining all this we obtain 298592 facets of HYPg in 86 
orbits. 

□ 

Definition 2. (i) Two Delaunay polytopes Pi and P2 are called isomorphic 
if there exists an isometry transforming one into the other. 

(ii) The automorphism group Aut(P) of a Delaunay polytope P is the set 
of all isometries leaving it invariant. 

Denote by Vi and V2 the vertex-sets of Pi and P2, respectively. In order 
to test the existence of such an isometry, it suffices to test if there exist a 
mapping <f> : Vi — > V2 satisfying to d((j>(x) , 4>{y)) = d(x,y) for all x,y £ Vi. 
The test of existence of an isometry is therefore a combinatorial problem. If the 
set of possible pairwise distances is {d\, . . . ,dh}, then the isomorphy problem 
becomes a problem of isomorphy of association schemes. The problem of 
computing the automorphism group of a Delaunay polytope is the same as of 
computing the automorphism group of an association scheme. Unfortunately, 
we do not know about a program doing such computations. 

Therefore, it will be presented here some algorithms doing computations 
using the nauty program (see (13)). Problem (ii) can be solved easily by 
associating, to every possible distance d between two vertices of P the graph 
Gd, formed by making two vertices adjacent if d(x,y) = d. So, we have 

Aut(P) = rf =1 G di , 

Problem (i) can be solved by associating to Pi and P2 a graph encoding pair- 
wise distances, but this approach leads to graphs, whose size is too big for 
the nauty program. So, we use procedures, which are not guaranteed to work 
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in any cases but which have worked so far in all cases considered of extreme 
Delaunay polytopes. 

Take two Delaunay polytopes Pi and P2 and compute their skeleton graphs 
(i.e. the graph formed by their vertices with two vertices adjacent if they 
generate a two dimensional face), G± and G2, respectively. If those graphs 
are not isomorphic, then P\ and P2 are not isomorphic. If they are isomor- 
phic, then there exists an isomorphism (f> : G\ — > Gi- If this mapping satisfies 
d((j)(x), 4>{y)) = d(x,y), then P\ and P2 are isomorphic. This approach does 
not always work, i.e. some automorphisms of G\ do not correspond to isome- 
tries of Pi and therefore the mapping (f> computed is not the right one. 

A similar approach is to consider the ridge graph (i.e. the graph formed by 
their facets with two facets being adjacent if their intersection is of dimension 
n — 2) G[ and G' 2 of Pi and P2. If G[ and G' 2 are not isomorphic, then 
Pi and P2 are not isomorphic. If they are isomorphic, then there exists an 
isomorphism (j)' : G\ — > G' 2 . Since every vertex is an intersection of some 
facets, this mapping lifts to an isomorphism <f> : G\ — > G2. We then test if <p 
defines an isometry. In all cases considered, (j) was an isometry, i.e. we were 
able to check the isomorphy of extreme Delaunay polytopes efficiently. There 
is no reason to think that this will always be the case. 



2. The Adjacency Method 

We consider in this section the details of the adjacency algorithm that takes 
an extreme ray e of HYP n+ ±, corresponding to an n-dimensional extreme 
Delaunay polytope and finds the adjacent extreme rays (ej)i<j< m in HYP n+ \. 

The problem in doing this computation is that the description of facets of 
HYP n is known only for n < 7 (see (9), (10) and (3)). Anyway, using the 
complete list of facets is not a good idea, since, for example, HYPg has at 
least 298, 592 facets. 

The algorithm, which we use, is as follows: 

(1) Given an initial extreme ray e, find the list of hypermetric vectors 
(pi)i<i<m such that H(bi)e = 0. We then take as initial list F of 
facets all vectors (6i)i<j< m plus all permutations of (l 2 , —1, n-2 ) (they 
correspond to triangular inequalities). 

(2) Define the cone C(F) by taking all hypermetric inequality in F. 

(3) Since e is an extreme ray of the cone C(F), we can find the extreme 
rays (ej)i<j< P of C(F) that are adjacent to e. 

(4) Test if ej is hypermetric, using the closest vector problem. 

(5) If the ray ej is hypermetric, then we finish. If not, then some hyper- 
metric inequalities H(b)d < are violated. So, we add those b to F 
and go back to step 2. 

Since the cone HYP n+ \ is polyhedral, the inner loop will eventually finish. 
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The key step in above algorithm is the ability to check if H(b)d < is true 
for all vector b <G Z n+1 with J2i h = 1- This is equivalent to 

n 

H(b)d B = 2(r 2 - Wvo-c + ^Tbilvi -v Q ) || 2 ) <0. 

i=i 

If we denote by L the lattice generated by the family (vi — t>o)i<i<n, then the 
problem is expressed as 

Does there exist v G L such that \\vq — c + v\\ < r? 

This problem is a closest vector problem (i) except that the scalar product 
(which comes from the distance d) is not always positive definite. 

We associate to d the Gram matrix G = (gij) defined by gij = \(dio + djo — 

If G admits a negative eigenvalue, then one can find a vector v G L, such 
that \\vo — c + v | < 0, thereby solving the problem. 

If G is positive but not positive definite, then this means that L is of dimen- 
sion inferior to n and that the n + 1 vectors form a Delaunay polytope of a 
lower dimensional lattice, which is not a simplex. In that case the algorithm 
is as follows: 

(1) Find the rank r of the matrix G. 

(2) Find a family F = {vi , . . . ,Vi r }, such that all points Vi can be ex- 
pressed in terms of F with integer coefficient (this is not always possi- 
ble). 

(3) If preceding step has succeeded, then the problem takes the form 

Does there exist hi-, such that || Y^j=Q^i 3 v ij ~ c \\ ^ r<!> 
The corresponding Gram matrix is positive definite. Therefore, this is 
a closest vector problem of type (i). 

3. Example of application of the Adjacency Method 

In (6) (see also (7)) an extreme 15-dimensional Delaunay polytope with 135 
vertices, i.e. having the minimal number of vertices, is given. We applied 
the adjacency method to one of its affine basis and found one new extreme 
Delaunay polytope, let us denote it by EDg of dimension 8. This polytope 
belongs to an infinite serie (ED n ) n >Q of extreme Delaunay polytopes (see 
(10)); this serie was found from EDg. 

We then apply the adjacency decomposition method to the polytope EDg 
and found 24 extreme Delaunay polytope of dimension 8. 

It is well-known (see (7)) that the Gosset polytope Gos is an extreme 7- 
dimensional Delaunay polytope. All orbits of affine basis of Gos were found 
by direct enumeration. 

We apply the Adjacency Method to all orbits of affine basis of Gos and 
found an extreme Delaunay polytope of dimension 7 with 35 vertices and a 
symmetry group of size 1440. In fact, this extreme Delaunay polytope was 
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already found by Erdahl and Rybnikov in (12). We conjecture that there is 
no other extreme Delaunay polytopes in dimension 7. 
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